Charge Density Waves and Axion Strings from Weyl Semimetals 
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We study dynamical instability and chiral symmetry breaking in three dimensional Weyl semimetals, which 
turns Weyl semimetals into "axion insulators". Charge density waves (CDW) is found to be the natural conse- 
quence of the chiral symmetry breaking. The Goldstone mode of this charge density wave state is identified as 
the axion, which couples to electromagnetic field in the topological 6E ■ B form. Our main finding is that the 
"axion strings" can be realized as the (screw or edge) dislocations in the charge density wave, which provides 
a simple physical picture for the elusive axion strings. These axion strings carry gapless chiral modes, there- 
fore they have important implications for dissipationless transport properties of Weyl semimetals with broken 
symmetry. 
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Introduction. Topological insulators are among the most 
active research fields in condensed matter physics recentlj^— . 
Among the remarkable aspects of topological insulators is the 
ubiquitous role played by Dirac fermions. In fact, most of 
the recently discovered topological insulators can be regarded 
as massive Dirac fermion systems with lattice regularization''. 
When the mass vanishes, we have massless Dirac fermions, 
which are two copies of Weyl fermion with opposite chirali- 
ties. The Dirac fermions, which obey the Dirac equation, are 
described by spinors with four components, while the Weyl 
fermions are two-component fermions described by the fol- 
lowing simple Weyl equation 



cr-y^ip- Elf/ 



(1) 



where + is referred as chirality (+ for left handed, - for right 
handed). Since the Dirac fermion can be decomposed into 
two copies of Weyl fermion, the latter is more elementary 
building block in the real world and theoretical models. In 
fact, it is our current understanding of nature that the elemen- 
tary fermions such as quarks and electrons fall into the Weyl 
fermion framework because the left-handed and right-handed 
fermions carry different gauge charges in the standard model 
of particle physics"^. 

It is worth noting that in writing down Eq.([T]i we have as- 
sumed that the two "Weyl points", at which the energy gap 
closes, are both located at k = 0. In the particle physics 
context, this assumption seems to be natural, however, in 
condensed matter physics, without imposing symmetry con- 
straints such as time reversal symmetry and inversion symme- 
try, Weyl points are generally located at different points in the 
momentum space. This fact has interesting consequences for 
Weyl fermion systems with broken symmetry, as we will show 
in this paper. 

Although Weyl fermion plays a crucial role in the de- 
scription of elementary fermions in nature, it has been stud- 
ied in the condensed matter context only very recently^^— . 
Weyl semimetals in three dimensions (3d) are analogous to 
graphene^^ in two dimensions (2d) in the sense that both are 
described in terms of gapless fermions with approximately 
linear dispersion, but the 3d Weyl semimetals are richer in that 
they are more closely related to various fundamental phenom- 
ena such as the chiral anomal y'"-"-^' . Unlike the topological 



insulators, whose transport is dominated by topologically pro- 
tected surface states, in the Weyl semimetal the bulk transport 
is most important. Their unique semi-metallic behaviors in 3d 
can potentially be engineered for semiconductor industry. 

Interaction effect plays a fundamental role in the dynam- 
ics of Weyl fermions. One possibility is the pairing inter- 
action which leads to the superconducting instability. Qi, 
Hughes and Zhang's fermi surface topological invariant^^ im- 
plies that if the pairing amplitudes have opposite sign for 
Weyl points with the opposite chirality, topological supercon- 
ductors are obtained. Another consequences of interaction, 
which we will focus in this paper, is the spontaneous chiral 
symmetry breaking. Chiral symmetry breaking is the phe- 
nomenon of spontaneous generation of an effective mass of 
Weyl fermions, namely a paring between the fermions (elec- 
trons) and anti-fermions (holes) with different chiralities. Due 
to chiral anomaly, the Goldstone mode 6 is coupled to electro- 
magnetic field as • B, therefore, this Goldstone mode is an 



axion 



'24-27 



In this paper, we focus on the chiral symmetry breaking in 
Weyl semimetal and axion strings in condensed matter con- 
text, which have new features absent in particles physics. We 
would like to mention that axion string can be realized on 
the surface of topological insulators with a magnetic domain 
wall^**, and axionic dynamics has been studied in topologi- 
cal magnetic insulators-'''^'^. In the present paper we study 
a new route to axionic dynamics through chiral symmetry 
breaking induced by interaction effect. The resultant states 
are charge density wave (CDW) states, which are experimen- 
tally observable. One of our main findings is that the (screw 
or edge) dislocations of CDW are exactly the "axion strings", 
which are important topological defects carrying gapless chi- 
ral modes. In these chiral modes electrons move solely in one 
direction without backscattering. In particles physics, axion 
strings have interesting cosmological implications such as the 
gravitational lenses effect^', but observable evidences of such 
axion strings are elusive. Axion strings in condensed matter 
systems have the advantage that they are much easier to de- 
tect. In the Weyl semimetal studied in the present paper, axion 
strings have important effects on the transport properties since 
they provide chiral modes supporting dissipatonless transport. 

Dynamical chiral symmetry breaking. We consider 
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the simplest model for dynamical chiral symmetry break- 
ing, which nevertheless captures the most salient physi- 
cal consequences. The Hamiltonian is given as H = 
H() + H\, in which //o is the free quadratic Hamil- 
tonian given as Hq = Zko-/; c^ff'^ff/sC^)'^'^' = 

-5l^k,k',q.?(q)(4+qa'^kff)(4'-q/j'^k'/?) = "5 Zq §(q)PqP-q. 

where a,p refer to all degrees of freedom other than the mo- 
mentum k, - '^k+qt'k, and -g(q) is the electron-electron 
interaction potential in the momentum space. 

Let us first turn off the interaction term Hi . Suppose that the 
free Hamiltonian Ho has two Weyl points located at Qi and 
Q2. We shall not be concerned with the full form of Hamil- 
tonian Hq since we mainly focus on the low energy physics, 
therefore we can expand 



h(k) = +Vf o- ■ (k - Q,) 



(2) 



where the + sign correspond to Q, with / - 1,2 respectively, 
cr - icTx, (Ty, CTj) are three Pauli matrices. Several model 
Hamiltonians for Hq can be found, for example, in Ref.^ ". 

For simplicity, in this paper we consider the case g(Q) > 0, 
where Q = Qi - Qa- Taking the mean field value of Hi, we 
can easily see that a chiral condensation v — 2k(4''''+Q) - 
Zk,a<4a'^k+Qa) leads to an energy gain -\g{Q)\v\^ < 0. This 
form of chiral condensation is a paring of electrons and holes 
with total momentum Q, which is reminiscent of the Larkin- 
Ovchinnikov-Fulde-Ferrell (LOFF) state of pairing between 
electrons with nonzero total momentum. More precisely, we 
can perform a Hubbard-Stratonovich transformation as 

* fpi<A 

where CL,p = CQ,+p and c^ p - CQ,+p, V is the volume of the 
system, and A is a momentum cutoff below which the linear 
Weyl fermion approximation is valid. Since we consider low 
energy dynamical symmetry breaking with m << A, contri- 
bution of fermions with |p| > A to the effective action of m is 
negligible. Such limit we take, which does not qualitatively 
affect the physical conclusions, will simplify our calculation 
significantly. Integrating out fermions, we obtain the effective 
action of m at the mean field level as 



' MF 



VP 



-TrlndetM 



(4) 



where V is the system volume, P - I /Temperature is the pe- 
riod along the imaginary time, which will be taken to be very 
large since we are concerned with the zero temperature limit. 
The fermion matrix M is given by 



M : 



St — Vfcr ■ p m 

m* dr + VfCr ■ p 



(5) 



It follows that the saddle point equation in the zero tempera- 
ture limit is 



1 r dojcPp 



1 



(xP- + v'|p^ + Imp 



(6) 



where g is the abbreviation of giQ). The solution to Eq.® is 
given as 



1 1 



1 



gc g Stt^v^ 



•In: 



A2 + Imp 



where 



.2,3 



8c 



(7) 



(8) 



Because we are concerned with the cases with \m\ « A, 
Eq.(|7]l is approximated by 



1 1 



1 



Imp In ■ 



A^ 



(9) 



It can be seen from Hq.® that the coupling constant must ex- 
ceed a critical value gc to give a physical meaningful solution 
for |m|. Therefore, the interaction should be sufliciently strong 
to lead to dynamical symmetry breaking (or called "exciton 
condensation"). 

Axion dynamics and topological theta term. We have seen 
in the previous section that when g > gc, chiral symmetry is 
spontaneously broken. The Ginzburg-Landau effective action 
of m can be expressed as 

5„, = J dtd\(^ydi,m*&'m + 6\mf + r]\mf) (10) 

Let us write m = |m| exp(iff). In the symmetry breaking phase, 
6 < and |m| develops an nonzero expectation value. From 
the Ginzburg-Landau effective action we can see that the fluc- 
tuations of |m| are gapped, while the 6 fluctuation is a gapless 
Goldstone mode. The Ginzburg-Landau action 5, „ is invariant 
under the transformation 6 9+ AO, but in fact this symmetry 
is broken by chiral anomaly, which endows a topological theta 
term to the effective action, as we explain as follows. We can 
perform a chiral transformation i^i — > e'^^-t//L, iffR — > e^'^^^i^^, 
then the transformation of the dynamically generated mass is 
readily obtained as m — > me^'^ = |m|. Although the action 
of low energy fermions is invariant under such chiral transfor- 
mation, the integral measure is not^^, thus the chiral transfor- 
mation generates an anomalous term asM 



* anomaly 



— r 

- r 



dtd^x6e''^^PF„yF 



fiv^ Ap 



dtd^xOE ■ B 



(11) 



where we have used the unit with h-c- \. 

Taking the above chiral anomaly into account, the gapless 
fluctuations of 9 is described by the following simple axionic 
effective action 

5s = y J" dtd^x(df,9f + J dtd\9E ■ B (12) 

where the notation "/a" is deliberately chosen because it is 
analogous to the pion decay constant f„, namely that f„ is the 
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"axion decay constant". We can also define a normalized field 
a - faO and put Eq.(fT2li into another form 

5„ = i J" dt£x(d^af + ^ J dtd\jE ■ B (13) 

The analogous effective action of the second term for pion is 
responsible for the famous two-photon decay of neutral pion. 
From the Ginzburg-Landau action in Eq.lfTOl) we can see that 
fa ~ \m\, therefore, we have the counterintuitive conclusion 
that when the chiral condensation goes to zero, the axion- 
photon coupling becomes strong. 

Various topological responses can be calculated from the 
effective action given in Eq.(fT2t. Taking derivative with re- 
spect to A^, we have the current 



(14) 



Such kind of responses have been studied in topological 
insulators^, and the case of Weyl semimetal is analogous in 
this respect. Let us consider the special case of a constant 
magnetic field B z along the z direction, then the charge den- 
sity given by Eq.(fT4l) is 



(15) 



This is consistent with the following physical picture. In a 
constant magnetic field Z?;Z, the motion of electrons in the 
xy plane is quantized into Landau levels with density of state 
eB./2n, then the problem is reduced to one dimensional, but 
we know that the charge density in one dimension can be 
found using the Goldstone-Wilczek formula^ - -^e^^dyO, 
therefore, the three dimensional charge density is obtained as 
/ = {eBJ2n){d£l2n) = ^,8,68,. 

Phase of charge density wave is the axion. Since we are 
concerned with low energy physics which is dominated by the 
two Weyl points, following the momentum space cutoff used 
in the previous sections, we expand the fermion operator as 
Cr - e'^^^cir + e'^-'cR r + ■ ■ where cr/i are cut off in the 
momentum space at A, i.e. c^/r - Zi|p|<A e''"^CL/fi.p + ■ ■ ■. and 
the "• ■ •" terms are high energy modes with |p| > A. The 
charge density is given by 



p(r) = clcr = 4rCL,r + c'^^Cr^t + (4,r'^«.>-« + ^.r^^.r^ 016) 

The mean field expectation value of chiral condensation is v - 
<c|^CR_r) - V"' lip(c|pCR,p) = m/g, which can be obtained 
from Eq.Q, therefore we have the charge density wave 



(p(r)) = po + vcos(Q-r-0) 

\m\ 

- po H cos(Q - r - 6) 



(17) 



where po is the constant background, while the cosine term 
indicates that the chiral symmetry breaking induces a charge 
density wave state with wave vector Q. 

When the CDW wavelength 2;t/|Q| is commensurate with 
the lattice constant, the chiral symmetry breaking transitions 



can be thought of as generalized Peierls transitions, which in- 
duces dimerization, trimerization, etc. Unlike the infinites- 
imal Peierls instability in one dimension, the generalized 
Peierls transitions in three dimensional Weyl fermion systems 
can occur only when the coupling is stronger than a critical 
value [see Eq.®]. The vanishing density of state at the Weyl 
point is the reason for the difference. 

The natural question is how to experimentally detect the 
CDW. It can be detected using bulk measurement, but it can 
also be detected by simpler surface measurements. The sur- 
faces naturally inherit CDW from the bulk, therefore the CDW 
can be detected using scanning tunneling microscope (STM). 
Denoting the angle between the surface normal and Q as a, we 
can obtain that the surface CDW wavelength as A2ci - jrJ^ — ,. 
When a 0, A2d oo, i.e. in the cases with Q perpendicular 
to the measured surface (a - 0), there is no surface CDW. 

It is worth noting that the interaction effect and CDW was 
studied in 2d Dirac systems in Ref.^^"^^. More recently, in- 
teraction effect on the surface of weak topological insulators 
has been studied in Ref.— , in which CDW has also important 
physical consequences. The relation to chiral anomaly is ab- 
sent in these studies because the systems considered there are 
in 2d, where the concept of chirality is lacking. 

Dislocations in charge density wave are axion strings. Let 
us turn to the central part of this paper, namely the identifica- 
tion of CDW dislocations as axion strings, which may provide 
dissipatonless chiral transport channel in 3d bulk materials. 
An axion string I is a one dimensional dislocation of axion 
field, around which the axion field changes by 27t, namely 
that 



X 



d0 = 27T 



(18) 



where C is a small contour enclosing / clockwise, is the ax- 
ion field coupled to the gauge fields in the form of ■ B. 
Axion strings are closely related to chiral anomaly, as was 
studied long ago in the work by Callan and Harvey^ in the 
particle physics context. In the Weyl semimetals studied in 
the present paper, the axion strings have clear geometrical pic- 
ture because is exactly the phase of CDW. More explicitly, 
suppose that Q = (Q^, Qy, Q,) = (0, 0, 0, then the maxima of 
charge density (peaks of CDW) is located at 2d planes (x, y, Zn) 
with 



+ 27m 

Q 



(19) 



where n ^integer When is shifted, the location of maximal 
density plane follows the shifting of 0. In fact, the shifting of a 
maximal density plane around the small loop C enclosing the 
axion string is readily obtained from Eq.(fT9]l as 



X 



dZn 



Q Q 



2n 



(20) 



which is exactly the wavelength of the CDW. The Burgers 
vector of the axion string as a dislocation of CDW is exactly 
(0, 0, 2nl Q). For a general CDW wave vector Q, the Burgers 
vector is readily obtained as 27rQ/|Qp. 



4 




(a) (b) 



FIG. 1: Axion strings as dislocations of charge density wave. 
(a)Screw dislocation. (b)Edge dislocation. The delineated sheets are 
the peaks of CDW. In both (a) and (b), the axion string is along the 
z direction. The Burgers vector is parallel to the axion string in (a), 
while perpendicular to it in (b). 



Let us refer the orientation of the axion string / as 1. Ac- 
cording to the relative orientation of 1 and Q, the axion string 
appears as different types of dislocation. When 1 is parallel 
with Q, we have screw dislocation [Fig. (a)], while when Q is 
parallel with 1, we have edge dislocation [Fig.(b)]. 

There are chiral modes propagating along the axion 
strings-'^, therefore, axion strings may serve as unique trans- 
port channels in a 3d materials with axionic dynamics. Such 
chiral modes carry dissipationless current just like the quan- 
tum Hall edge and quantum anomalous Hall edge states, but 
the former are more robust because they are buried in 3d bulk 
materials rather than placed at 2d sheets. 

It is worth mentioning that dislocations in topological ma- 
terials has also been studied in Ref.^^, but we would like to 
emphasize several prominent differences between the axion 
strings studied here and the dislocation lines in weak topolog- 
ical insulators studied in Ref.^^. First, in Refi^ the disloca- 
tions carrying gapless modes are indeed dislocations of crys- 
tal lattice, while in our paper the crystal lattice remains intact, 
and axion strings are just dislocations of CDW. Secondly, the 
gapless modes studied in Ref.^*^ are helical modes, which are 
unstable towards back scattering if time reversal symmetry is 
broken, while the gapless modes living on the axion strings 
studied in the present paper are robust chiral modes. It is also 
worth noting that in Ref. ' ' line dislocation with chiral modes 
was studied, but CDW is absent there, more importantly, the 
bulk is also gapless there and the coupling between disloca- 
tion mode and bulk mode can induce dissipation. 

To conclude this section we remark that the formation of 
axion strings in Weyl semimetals can be triggered by rapidly 
lowering the temperature from T > Tc to T < Tc, where is 
the critical temperature of chiral condensation (Kibble-Zurek 



mechanism). 

Discussions and Conclusions. It is worth mentioning that 
an attractive interaction -giQ) < at momentum Q is needed 
for the chiral symmetry breaking. The values of giQ) for var- 
ious materials depend the material details, which is not within 
the scope of this paper. It is useful to mention that an ef- 
fective attractive electron-electron interaction can appear at 
some special momenta commensurate with the reciprocal lat- 
tice. Such electron-lattice coupling effect is responsible for 
the Peierls transitions in one dimensional systems, therefore 
we see that a sufficiently strong electron-lattice coupling is 
able to cause dynamical chiral symmetry breaking if Q is 
commensurate with the reciprocal lattice. 

Our model provides a geometrical picture of axion, which 
manifests itself as CDW. CDW phases have (edge and screw) 
dislocations, which are exactly the axion strings. The CDW 
and associated dislocations are much more visible than the 
elusive axion and axion strings in previous settings. One of 
our main results in this paper is the identification of axion 
strings as CDW dislocations, which has no analog in particle 
physics. The axion strings have one dimensional robust chiral 
modes along them, which have great potential applications if 
the chiral symmetry breaking (exciton condensation) of Weyl 
fermion is realized in experiment. In this paper we studied the 
general cases with Q ?t 0. When Q = there is no CDW 
associated with the chiral symmetry breaking, but the axion 
strings do exist and have important implication for 3d trans- 
port properties. 

To conclude this paper, let us mention that the dynamical 
chiral symmetry breaking in Weyl semimetals discussed in 
this paper can also be realized in cold atom experiment. Al- 
though the anomalous 0E ■ B coupling is absent because atoms 
are neutral, the atom density wave and its dislocations can be 
realized. Because the flexible tunability in cold atom exper- 
iments, it is interesting to realize Weyl fermions and chiral 
symmetry breaking in cold atom systems. 
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Note added. After this work was finished, we became aware 
of a related work**" on symmetry breaking in Weyl semimetal 
by Zyuzin and Burkov, though CDW and axion strings were 
not studied. Due to nonzero density of states considered in 
their work, the gap equation is different from ours and an in- 
finitesimal interaction will induce chiral condensation in that 
case. 
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